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Abstract 

In this paper, we present an interference model for cognitive radio (CR) networks employing power 
control, contention control or hybrid power/contention control schemes. For the first case, a power 
control scheme is proposed to govern the transmission power of a CR node. For the second one, a 
contention control scheme at the media access control (MAC) layer, based on carrier sense multiple 
access with collision avoidance (CSMA/CA), is proposed to coordinate the operation of CR nodes 
with transmission requests. The probability density functions of the interference received at a primary 
receiver from a CR network are first derived numerically for these two cases. For the hybrid case, 
where power and contention controls are jointly adopted by a CR node to govern its transmission, 
the interference is analyzed and compared with that of the first two schemes by simulations. Then, 
the interference distributions under the first two control schemes are fitted by log-normal distributions 
with greatly reduced complexity. Moreover, the effect of a hidden primary receiver on the interference 
experienced at the receiver is investigated. It is demonstrated that both power and contention controls 
are effective approaches to alleviate the interference caused by CR networks. Some in-depth analysis 
of the impact of key parameters on the interference of CR networks is given via numerical studies as well. 

Index Terms - Cognitive radio, interference modeling, hidden primary receiver. 



IEEE TRANSACTION ON COMMUNICATIONS, VOL. XX, NO. Y, MONTH 2010 



2 



I. Introduction 

With the requirement to improve spectrum utilization, the newly emerging cognitive radio 
(CR) technology [Q]|- Hi has attracted increasing attention. A CR network is envisioned to be 
capable of reusing the unused or underutilized spectra of incumbent systems (also known as pri- 
mary networks) by sensing its surrounding environment and adapting its operational parameters 
autonomously. A CR system may coexist with a primary network on either an interference-free 
or interference-tolerant basis [5J. For the former case, the CR system only exploits the unused 
spectra of the primary network, which consequently guarantees no interference to primary users. 
For the latter case, the CR system is allowed to share the spectra assigned to the primary network, 
under the condition that the CR network must not impose detrimental interference on the primary 
network. Therefore, modeling and analyzing the interference caused by CR networks is of great 
importance to reveal how the service of a primary network is deteriorated and how CR networks 
may be deployed. 

In the literature, the existing research on interference modeling for CR networks mainly 
falls into three categories: spatial, frequency-domain and accumulated interference modeling. 
For spatial interference modeling, the fraction of white spaces available for CR networks was 
investigated in [|6l and 0. In 0, the region of interference for CR receivers and region of 
communication for CR transmitters were studied for the case where a CR network coexists 
with a cellular network. The interference from CR devices to wireless microphones operating 
in TV bands was analyzed in J9]|, where the loss of reliable communication area of a wireless 
microphone due to the existence of CR devices was examined. CR interference in the frequency 
domain was also researched in the literature, e.g., the interference due to out-of-band emission 
of a wireless regional area network (WRAN) was analyzed in ifTOl . 

As for accumulated interference modeling, in 0TJ, the aggregate interference power from a 
sea of CR transmitters surrounding a primary receiver was derived. Also, the accumulated CR 
transmission power perceived at a primary receiver was given by integrating over the "CR sea" 
with a certain power density. The performance of a primary system was evaluated in [[T2| in terms 
of outage probability caused by the interference from CR networks. The outage probability was 
derived for both underlay and overlay spectrum sharing cases. In [fT3l the aggregate interference 
from multiple CR transmitters following a Poisson point process was approximated by a Gamma 
distribution and the probability of interference at a primary receiver was also given. It is worth 
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noting that only pathloss was assumed for the interfering channel in IflTTl - |fT3l . Their work was 
extended by taking both shadowing and fading into account in lfl4ll and [fT51 . Moreover, the 
probability density function (PDF) for accumulated interference and outage probability due to 
the aggregate interference from CR nodes were also derived in |[T4ll and Ifl51 , respectively. 

However, in all the previous works (6]|- [fT51 . the CR transmitters were assumed to transmit 
at a fixed power level, i.e., no power control for CR transmitters was considered. Moreover, 
the CR nodes were all assumed to communicate with each other simultaneously. Thus, no 
contention control scheme was employed at the cognitive media access control (MAC) layer. 
Some preliminary results on CR interference modeling were obtained in [fT6ll by incorporating 
either power or contention control scheme. In this paper, we extend the aggregate interference 
modeling in the following aspects. Firstly, a more realistic power control scheme than that in lfl6l 
is proposed, and a hybrid power/contention control scheme is introduced. Secondly, the PDFs 
of interference perceived at a primary network from a CR network are derived numerically for 
the cases of power or contention control. The interference distribution of the hybrid control 
scheme is also analyzed and compared with that of the pure power control and pure contention 
control schemes by simulations. Furthermore, for the power and contention control schemes, 
their interference distributions are fitted by log-normal distributions, which greatly reduces 
computational complexity compared to a numerical approach to obtain PDFs. Finally, the impact 
of a hidden primary receiver on the aggregate interference is investigated for all the three schemes. 
The impact of several key parameters on the resulting interference is evaluated as well, which 
provides some insights for the deployment of CR networks. 

The remainder of this paper is organized as follows. The system model is elaborated in 
Section II. The detailed interference modeling is presented in Section III. In Section IV, the in- 
terference distributions are approximated by log-normal distributions. We incorporate the hidden 
primary receiver problem in Section V. The impact of several key parameters on the interference 
is analyzed via numerical studies in Section VI. Finally, Section VII concludes the paper. 

II. System Model 

The system model is illustrated in Fig. 1. It consists of a CR network coexisting with a 
primary transmitter-receiver pair. The active CR transmitters are distributed in a 2-dimensional 
plane outside the interference region (IR) of the primary receiver as shown in Fig. 1. The IR 
is a disk centered at the primary receiver with a radius R. CR transmission is forbidden within 
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this circular IR in order to protect the primary receiver against co-channel interference from the 
surrounding CR transmitters, since it is assumed that all the CR transmitters reside in the same 
frequency spectrum as the primary transmitter. We model the aggregate interference received at 
the primary receiver due to the existence of a CR network and investigate the impact of CR 
network deployment parameters on the resulting aggregate interference. 

The underlying interference channels from CR transmitters to the primary receiver experience 
pathloss, shadowing and fading. The pathloss function g(rj) is 

g ( rj ) = rj 13 (1) 

where rj is the distance between the jth (j = 1, 2, • • ■ ) active CR transmitter and the primary 
receiver and (3 is the pathloss exponent. The composite model for shadowing and fading can 
be expressed as the product of the long term shadowing and the short term multipath fading. 
In this paper, log-normal shadowing and Nakagami fading are considered. Let hj denote the 
channel gain for the composite shadowing and fading of the interference channel from the jth 
active CR transmitter to the primary receiver. The PDF of the composite channel gain hj can 
be approximated by the following log-normal distribution ifTTl 



where the mean fi and variance a can be expressed as 

k 



m— 1 j \ 

£ --ln(m)- 0.5772 +^ (3) 



° 2 = Y, , 1 n 2 +^ W 

with m standing for the Nakagami shape factor and fin and a n denoting the standard mean and 
variance of the log-normal distribution, respectively. 

Let pj denote the transmission power of the jth active CR transmitter. The accumulated power 
of the instantaneous interference received at the primary receiver can be expressed as 

oo 

y = 12Pi9(rj)h s . (5) 

In this paper, we investigate the characteristics of the aggregate interference from all CR transmit- 
ters employing the following three different schemes: (i) power control, (ii) contention control, 
and (iii) hybrid power/contention control. 
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A. Power Control 

In this scenario, the distribution of active CR transmitters follows a Poisson point process 
with a density parameter A for the density of CR transmitters on the plane. 

The transmission power of a CR transmitter is governed by the following power control law 

/ \ ( r ~ ) -Pmaxj < r cc . < T pwc 

P P wc(r C c 3 .) = < V pvcJ (6) 

where r CCj is the distance from the jth active CR transmitter to its nearest neighbouring active 
CR transmitter, a is the power control exponent, P max is the maximum transmission power for 
CR transmitters, and r pwc is the power control range, which determines the minimum r CCj leading 
to maximum CR transmission power P max - Compared to the power control law in lfT6l . a new 
parameter r pwc is introduced here to adjust the range of the power control. We assume that 
the power control exponent a is equal to the pathloss exponent in (OQ) throughout the paper. 
The above proposed power control scheme is designed in such a manner that the interference 
caused by the jth active CR transmitter to its nearest active CR transmitter due to pathloss is 
P P wc( r cc :) )fi , ( r cc ; ,)- h is clear that within the power control range r pwc , this interference is equal 
to a constant P max /r^ wc . But beyond the power control range, the interference is less than that 
constant. In other words, at any CR transmitter the interference from the nearest neighbouring 
CR transmitter is capped and independent of the nearest neighbour distance within the power 
control range. It is worth noting that for each CR transmitter the information of its nearest 
neighbour distance is indispensable to determine its transmission power. Therefore, to facilitate 
the abovementioned power control scheme, either a central console having the global position 
information of all active CR transmitters or a distributed sensing scheme for CR transmitters like 
pilot sensing lfT9l is required, the detail of which is, however, beyond the scope of this work. 
When CR transmitters follow a Poisson point distribution with a density A, the PDF of r CCj can 
be given as ll20ll 

f cc (x) = 2nXxe~ Xnx2 . (7) 

B. Contention Control 

Unlike the previous power control scheme, for the case of contention control every active 
CR transmitter has fixed transmission power p, but their transmission is governed by contention 
control to determine which CR transmitters can transmit at a given time. We assume that the 
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multiple access protocol carrier sense multiple access with collision avoidance (CSMA/CA) 
is employed, like in IEEE 802.11 networks. Every CR transmitter senses the medium before 
transmission. If the medium is busy, namely, the CR transmitter detects transmission from 
other CR transmitters within its contention region, it defers its transmission. Otherwise, the 
CR transmitter starts its transmission. As a result of the contention control shown, all the active 
CR transmitters are separated from each other by at least the contention distance, which is the 
minimum distance rf min between two concurrent CR transmitters. 

The distribution of the active CR transmitters under the contention control can be modeled 
as a Matern-hardcore (MH) point process [fT8ll . which can be considered as a thinned process 
from a Poisson point process EUl . The thinning operation deletes some points from the original 
Poisson process under certain criteria. The MH process $ m h is the result of dependent thinning 
from a Poisson point process $, i.e., deleting or retaining a point depends on previous deletion 
operations. The mathematical expression of the MH process is given by ll20l 

&mh={x e$:m(x)<m(y) for all y in $ n C(x, d min )}. (8) 

Each point x in the original Poisson point process $ is marked with a random variable m(x) 
uniformly distributed in (0,1), while C(x, d m ; n ) is a disk centered at point x with the radius d m i n . 
The retaining probability q m h for the MH process, which is the probability of a point from a 
Poisson point process with a density A surviving the thinning process, is given by ll20l 

i ^,—Xird 2 ■ 

_L — C mm 

A7r «min 

C. Hybrid Power/Contention Control 

The aforementioned power control scheme regulates the transmission power of each CR 
transmitter according to its nearest neighbouring transmitter distance, while the contention control 
determines which CR transmitter can transmit at a time instant with fixed transmission power. 
A natural extension of the above two interference management schemes is to implement both 
schemes in the same system. This is termed hybrid power/contention control and it works in 
the following manner. The contention control scheme is first applied, resulting in a set of active 
CR transmitters following an MH point process. Then, a power control scheme similar to © is 
employed to adjust the transmission power of each active CR transmitter according to the distance 
to the nearest neighbouring active transmitter. The following power control law is adopted in 
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the hybrid control scheme 



Ph y b(r) 



1 fe) a * 



'mm 



< r < r hyh 



r > r hyb 



(10) 



where r is the distance from an active CR transmitter to its nearest neighbouring active CR 
transmitter, a is the power control exponent as in ©, and r hyh is the power control range 
similar to r pwc in © except that it also determines the maximum transmission power, i.e., 



yT^j P- I* * s °b v i° us that a larger r hyb leads to a larger maximum CR transmission power 
and, consequently, longer communication range for CR transmitters. The above power control 
law (flOl) guarantees that when a pathloss channel is considered for each active CR transmitter, the 
perceived interference caused by its nearest neighbouring CR transmitter is Phyb(r)g(r), which 
is (i) a constant p/d^ in within the power control range r hyb and (ii) less than the constant p/<i min 
when the distance r is larger than the power control range. 



We intend to model the aggregate interference from CR transmitters employing the three dif- 
ferent interference management schemes introduced in Section II by finding their corresponding 
PDFs. We apply the methodology used, for example, in lfl4ll and ETI to derive the PDFs. First, 
the characteristic functions of the interference under different system models are derived. Then, 
the PDFs of the aggregate interference are obtained by performing an inverse Fourier transform 
on their characteristic functions. 

A. Power Control 

When all the CR transmitters follow a Poisson point process distribution and employ the 
power control scheme proposed in (J6]), we can adopt the characteristic function-based method 
as in ffT4l . ETTl - Il23l and obtain the following characteristic function (py(co) of the aggregate 
interference Y at a primary receiver from all CR transmitters 



where / p (-) is the PDF of the transmission power p pwc (r CCj ) of a CR transmitter defined in © 




III. Interference Modeling 




(11) 



and 




(12) 
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In (fl"2l) . denotes the inverse function of g(-) in (OQ). For the derivation of (flTT) . the following 

fact is used: the distances from the jth CR transmitter to the primary receiver rj (J = 1, 2, • • ■ ) 
have independent and identical uniform distributions for a given number of CR transmitters [I2TTI . 
Their PDFs have the following form [12111 

2x/(l 2 -R 2 ), R<x<l 

~ ~ (13) 
0, otherwise 

when CR transmitters are distributed within an annular ring with inner radius R and outer radius 

/. In (fTT|) . p is a function of r cc as shown in ©, so the expectation of T{ujph) over p equals 

that of T(uip pwc (r cc )h) over r cc . Using the PDF of r cc given in CD), (fTTI) can be rewritten as 

0yO) = exp ( A7r / f h (h) f cc (r)T(up pwc (r cc )h)drdh) . (14) 

Moreover, (fl4"l) can be written as (see Appendix A for the detailed derivation procedure) 

iuir a Pmaxg{R)h 



<py(u)=exv{\>K I f h (h)J f cc (r) 



iur a P max h f 9{R) 2 i^ r "iw<. , 
+ / t T e r *™ a dt 



r 



pwc 







drdh 



r r°° rg(R) 2 

+\n f h (h) f cc (r) R\l-e* J °W p -~ h )+iujP^h t~?e M - h dt 

J H J r DWC J 



drdh, 



(15) 



Finally, we obtain the PDF of the interference by performing the inverse Fourier transform 

on 4>y{oj) as 

1 



My) 



2tt 



(f)Y(u)e- 2 ™ y du. 



(16) 



Equations ( fT5T ) and (fT6b serve as general expressions for the characteristic function and PDF, 
respectively, of the interference under the power control scheme. As a special case, when the 
pathloss exponent ft = 4 and the radius of the interference region R = 0, the PDF fy{y) can be 
further simplified through similar steps to that used in [21 J and obtained as 



7T 



f Y {y) = -K\y^' 2 e W 



ti 3 \ 2 K 2 



Mj 



where 



K=^KZ / f h (h)Vhdh 
Jh 



27crXe 



-Xnr 2 



dr + e 



-Atttv 



'0 V pwc , 

The detailed derivation procedure for K can be found in Appendix B. 



(17) 



(18) 
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B. Contention Control 

As mentioned in Section II.B, the distribution of CR transmitters can be modeled as an MH 
point process when the contention control is adopted. The MH process is a dependent thinning 
process from the original Poisson point process, which means that the positions of CR transmitters 
are correlated to each other. However, it is very difficult to obtain the distribution function like 
(fl"3l) for an MH point process in order to model the distance from an active CR transmitter to 
the primary receiver. Instead, we approximate the MH point process as an independent thinned 
Poisson point process with retaining probability q m h given by ©. Then, the transmission power 
for the jth CR transmitter is pj = {0,p}, which is a random variable taking values p or with 
probabilities p m ^ and 1 — p m h, respectively. To this end, the contention control scheme can be 
interpreted as follows: all the CR transmitters still follow the original Poisson point process with 
intensity A, but the jth CR transmitter has probability g m h to transmit at power level p. The 
characteristic function of the accumulated interference can be found as 



The detailed derivation of (fl~9l) is presented in Appendix C. 

Moreover, the PDF of the interference can be obtained from (fl9l) and (fT6l) . As a special case, 
when no IR is implemented and the pathloss exponent f3 = 4, this PDF can be simplified as 
(fTTT) with 



It is worth noting that the approximation for the MH point process actually ignores the 
dependence among the CR transmitters and treats an MH point process as a result of independent 
thinning process from an original Poisson point process. The accuracy of this approximation is 
evaluated in Section IV. 

C. Hybrid Power/Contention Control 

So far, the PDFs of the interferences received at a primary receiver from a CR network 
employing power control and contention control schemes have been derived. In order to model 
the aggregate interference under the hybrid control scheme, the nearest neighbouring distance 
distribution function analogous to © for active CR transmitters is indispensable to evaluate the 
transmission power designated in (flOl) . Unfortunately, there is no closed-form expression for 




(19) 




(20) 
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the nearest neighbour distance distribution function for an MH point process E4ll . Alternatively, 
several estimators have been used to statistically estimate the nearest neighbour distance distri- 
bution function in practice 11251 . However, statistical estimation is not practical for deriving the 
characteristic function in our case. Thus, we approach this problem numerically. 

The PDF for the aggregate interference under the hybrid control scheme is simulated in 
Fig. 2, where the interference PDFs for power and contention control are given as well for 
the purpose of comparison. It can be seen from this figure that both the mean and variance of 
the aggregate interference increase for the hybrid control scheme compared to either power or 
contention control schemes. However, the boosted interference is paid off by the increased CR 
communication area (coverage) for the hybrid control scheme. We define the coverage of each 
CR transmitter as a circular disk centered at a CR transmitter with radii being min(r/2, r pwc /2), 
d min /2 and min(r/2, r hyb /2) for power control, contention control and hybrid power/contention 
control schemes, respectively. Then, the received signal power at cell edge of a CR transmitter 
due to pathloss is 2 /3 P max /r pwc /3 , 2 l ^p/d^ ain and 2 /3 p/d^ lin for the above three aforementioned 
schemes. For the sake of comparison, let r pwc = d m - m and P max = p, which guarantees that the 
strength of the received signal power at cell edge of a CR transmitter is the same for all the 
three schemes. The overall coverage of the CR netwrok under different control schemes can 
be investigated numerically. With this setup, the overall coverage ratio for the power control, 
contention control and hybrid power/contention control is 1.0093, 1, and 2.0229, respectively. 
Two interesting facts are unveiled from this experiment. Firstly, the power control scheme 
leads to slightly smaller interference and slightly lager coverage compared to the contention 
control scheme, which suggests that power control is preferable to contention control in terms 
of lower resulting interference and larger coverage if the CR system can afford the complexity 
introduced by implementing the power control scheme. Secondly, the hybrid scheme tends to 
cause higher interference, but it greatly enlarges the coverage compared to power and contention 
control schemes. 

IV. Analytical Approximation 

In the previous section, to derive the PDFs for aggregate interference, the characteristic 
function-based method has been used which consists of two steps. Namely, characteristic func- 
tion computation and Fourier transformation. This interference modeling approach is extremely 
computation-intensive, since generally closed-form expressions are not admitted for either step 
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and the computations in both steps have to be performed numerically. It is desirable to model 
the aggregate interference with less complexity. An alternative approach to model the interfer- 
ence, which greatly reduces complexity, is to approximate interference PDFs as certain known 
distributions. Observations from Fig. 2 suggest that the interference distribution for either power 
or contention control is positively skewed and heavy-tailed, which suggests a log-normal dis- 
tribution. Thus, in this section, we fit the aggregate interference under power and contention 
control schemes to log-normal distributions. The theory behind the log-normal fitting is based 
on the following two facts. It has been shown that the sum of interference from uniformly 
distributed interferers in a circular area is asymptotically log-normal lfT5l . Il2~6ll . This ensures that 
the aggregate interference in these two schemes can be approximated as log-normal distributed. 
Meanwhile, the sum of randomly weighted log-normal variables can be modelled as a log- 
normal distribution as well 1127 1. which guarantees that the aggregate interference is still log- 
normal distributed even if the effect of shadow fading (O is taken into account. In what follows, 
the log-normal fitting is performed using a cumulant-matching approach ll28l . where the first 
two order cumulants of the aggregate interference Y in © are used to estimate the mean and 
variance of the log-normal distribution function. Therefore, the exact PDFs of interference can be 
obtained. Fortunately, these cumulants have closed-form expressions for both control schemes. 
Consequently, it significantly reduces the complexity compared to the interference modeling 
carried out in Section III. 

For the PDF of a log-normal variable x 

1 /-(\n(x) - u) 2 \ 

p(x) = ^=^exp 1 \[ 2 N (21) 



its mean \i and variance a 2 can be estimated using its first two order cumulants k\ and k 2 as 
follows [|29]|: 

fx = In kl (22) 



k 



a 2 = ln(j| + l). (23) 

In the context of interference distribution fitting, the nth cumulant k n of the aggregate interference 
Y can be obtained from its characteristic function (j>y(u) via the following equation 



kn — T- 



d n \n(j) Y (u) 



(24) 

J oj=0 
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A. Power Control 

From CES) and (|24l) . the cumulants for aggregate interference under the power control scheme 
can be derived as (see Appendix D for detailed derivation) 



(n/3 - 2)R n P~ 2 



na(na — 2) ■ • • 2 
r pwc ™(27rA)^ 

^na(na — 2 

i=l 



^ _ e -A7rr pwc ' 



net 



2i + 2) 



(27r\r pwc 2 y 



na-2i -\Trr pv ,c 
' pwc C 



(25) 



To evaluate the accuracy of the log-normal approximation for the power control case, some 
comparisons are performed in Fig. 3(a). It can be seen from Fig. 3(a) that there is fairly 
good agreement between the interference PDFs derived in Section III and the approximated 
counterparts. This approximation approach can be applied to both the pathloss-only and shadow 
fading channels. 



B. Contention Control 

Following the similar steps as in Appendix D and given the characteristic function (fT9l) for 
the aggregate interference under contention control and also using (T24l) . we can find the nth 
cumulant k n of aggregate interference as 



k 



Xnq, 



mh 



n 



-R 2 (ipg(R)h) n + n(iph)' 



g(R) 



-n-l-4 



dh 




g n -HR)-R 2 g n (R)\p n / f h (h)h n dh 



1 _ e - A7ra! min ) e n ^+ L 



(26) 



- 2)d£ ln iP*-* 

The accuracy evaluation of log-normal approximation under the contention control scheme 
is also performed and shown in Fig. 3(b). It can be seen from this figure that the log-normal 
approximation is fairly accurate compared to the simulated interference PDFs for either pathloss- 
only or shadow fading channels. Moreover, the derived interference PDF obtained from (fT6l) and 
(TT9l) is validated against the simulated counterpart in Fig. 3(b) as well, which suggests that the 
approximation for the MH point process in the derivation is reasonable. 
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V. Imperfect Primary System Knowledge 



In practice, some information about the primary system may not be perfectly known. One 
prominent example is the location of the primary receivers, which is usually required by CR 
networks in order to protect primary receivers from interfering CR transmitters. However, this 
information is not always available, especially in the case of passive primary receivers, i.e., when 
the primary receivers are hidden from CR networks. It is widely accepted that passive receiver 
detection techniques can be used or developed in the context of CR networks. For example, one 
of such primary receiver detection techniques is reported in OTTl . Nevertheless, its applicability 
is still not convincingly viable since it requires deploying sensor nodes close to primary receivers 
and much coordination is involved between these sensors and CR networks as well. The most 
commonly used and also the simplest approach to protect the primary receiver is to regulate the 
transmission of the CR network based on primary transmitter sensing, assuming that primary 
receivers are in close proximity to the primary transmitter. In this section, we evaluate the effect 
of a hidden primary receiver on the resulting interference to primary receivers. 

Consider a primary and CR coexisting systems depicted in Fig. 4, where an IR with radius 
R centered at the primary transmitter is introduced. All CR transmitters are distributed in the 
shaded concentric ring with inner radius R and outer radius I. Let 9 be the angle between 
the line joining the primary receiver and a CR transmitter and the line joining the primary 
transmitter-receiver pair. The distance from the CR transmitter to the primary transmitter is r 
and the distance between the primary transmitter-receiver pair is r p . Then, the distance between 
the CR transmitter and the primary receiver r cp can be expressed as 



where r is distributed as in (fT3l) and 9 is uniformly distributed in [0, 2tt] if a Poisson point 
process is assumed for the CR transmitter distribution. 

A. Power Control 

Under the power control scheme proposed in Section II.A and the system model given in 
Fig. 4, the characteristic function of aggregate interference 4>y(ui) can be written as follows (see 




(27) 
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Appendix E for the detailed derivation): 



V(w) = limexp(A f f h {h) I ^ f cc {x) [ [ e ioJ te) Pm ^ {x)9{rcp{r ' 9))h r- r dr d6 dx dh 
I Jh Jo ' Jo Jr 

p POO p2ll pi j 

+ ^ fh{h) fcc(x) e ^P^(x)g( rcp (r,e))h r _ r drdOdxdh} . (28) 



H J r pwc JO J R 



Applying the log-normal approximation method used in Section IV, we obtain the A;th cumulant 
of interference as 



k n = lim A (/ f h (h) f ^ f cc (x) [ * [ ( rap ™^ 9 M r > ^^ rdr d6 dx dh 

r poo r2n pi j 

+ h(h) f cc (x) / [P max (x)g(r cp (r,6))h} n rdrd6dxdh\. (29) 
Jh Jr pvlc Jo Jr J 

As can be seen from (|29l , unlike (1231) , the kth cumulant does not have a closed-form expression. 
However, the complexity of obtaining the exact interference PDF from (|29l is still smaller than 
that of the numerical method in Section III. 

An experiment is performed in Fig. 5(a) to examine the effect of hidden primary receiver 
on the resulting interference compared to the interference for the case of perfect knowledge of 
primary receiver location. We consider a pathloss-only channel in this figure. It can be seen from 
the figure that the hidden primary receiver problem boosts the interference in terms of increased 
interference mean and variance. This figure also shows that the log-normal approximation still 
fits well the interference distribution in this scenario. 

B. Contention Control 

Under the contention control scheme proposed in Section II.B and the system model given in 
Fig. 4, the characteristic function of aggregate interference 0y(o;) can be expressed as 

<Mw) = I™ exp {g mh AvrA (E (^™^ h ) - l) } 

I— >OQ 

= lim exp |g mh AvrA (J f b (h)J -|- J exp [iupg (r cp (r, 6))h] ^-dr dddh-\ 



= lim exp < g m hA / fh(h) / exp [iupg(r cp (r, 8))h] r — rdr d9 dh > , (30) 
'- >DO I Jh Jo Jr J 

with Di = l 2 -R 2 . 

Using the same log-normal approximation method as in Section IV, the A;th cumulant of 
interference can be written as 

p p2n pi 

k n = lim g mh A / f h (h) / / \pg(r cp (r,6))h] n r - rdr d6 dh. (31) 
Jh Jo Jr 
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The effect of hidden primary receiver under contention control is evaluated in Fig. 5(b), where 
a pathloss-only channel is assumed. As we can see from this figure, the uncertainty about the 
primary receiver location leads to interference with larger mean and variance as compared to that 
in the case with perfect knowledge of primary receiver location. Moreover, it can be seen from 
this figure that the log-normal fitting for the interference is fairly accurate and the approximation 
approach is still applicable in this scenario. 

For the case of hybrid power/contention control, the effect of hidden priamry receiver cannot 
be examined analytically because the closed-form interference PDF is not available. Therefore, 
it is analyzed numerically in Fig. 6, whose initial setup is the same as the one used in Fig. 5(b) 
except that the power control range is r hyh = 30 m. It can bee seen from Fig. 6 that the uncertainty 
about the primary receiver location boosts the interference in terms of increased mean, variance, 
and heavier tails for the hybrid control scheme as well. More interestingly, another two facts can 
be found by comparing Figs. 5 and 6: (i) the hidden primary receiver phenomonon has similar 
impact on the pure power and pure contention control schemes; (ii) the hybrid power/contention 
control scheme is less sensible to the phenomenon of hidden primary receiver than any of the 
other two schemes. 

VI. Numerical Studies 

The aggregate interference power from CR transmitters employing power control or contention 
control is investigated numerically in this section. For the power control scheme, Fig. 7(a) shows 
the effect of different power control parameters on their resulting aggregate interference. The 
detailed setup for the initial power control scheme is as follows: the maximum transmission 
power for each CR transmitter P max = 1 W, the density of CR transmitter A = 3 user/10 4 m 2 , 
the IR radius R = 100 m, the power control range r pwc = 20 m, the pathloss exponent (5 = 4 
and the power control exponent a — 4. From the two rightmost PDFs in this figure, it can 
seen that introducing power control scheme actually shifts the interference distribution leftwards 
compared to the distribution without power control. It means that the power control scheme 
can reduce the interference experienced at the primary receiver in terms of reducing its mean 
and slightly decreasing its variance. When deploying a CR network under the power control 
scheme, its resulting interference can be controlled by manipulating several parameters including 
-fmax> r P w« K and R. It can be seen in Fig. 7(a) that the interference can be reduced by either 
decreasing the maximum transmission power and/or CR density, or increasing the power control 
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range and/or IR radius. Interestingly, it also suggests that adjusting the IR radius is an effective 
way to control the interference, since the interference is more sensitive to the IR radius than to 
any other parameter as demonstrated in Fig. 7(a). Meanwhile, the interference is least sensitive 
to the CR user density in the sense that halving A leads to higher interference compared to 
doubling r pwc , halving P max or doubling R. 

For the contention control scheme, the impact of contention control parameters on the resulting 
interference is depicted in Fig. 7(b), whose initial setup is the same as that of Fig. 7(a) except that 
the transmission power for each CR transmitter is p = 1 W and the contention control range is 
^min = 20 m. It can be seen from the two rightmost PDFs in Fig. 7(b) that the contention control 
scheme results in an interference distribution with reduced mean like the power control scheme 
in Fig. 7(a). Meanwhile, the interference can be reduced by decreasing p, A, and/or increasing 
R or d min . It can be observed by comparing Fig. 7(b) with Fig. 7(a) that (i) increasing the IR 
radius is an effective approach to reduce the interference for both the power and contention 
control schemes. However, the power control scheme is more sensitive to the IR radius than 
the contention control one; (ii) reducing the transmission power and/or CR transmitter density 
affects the interference in the very similar manner for these two control schemes. 

Finally, the impact of shadow fading on the aggregate interference is investigated for dif- 
ferent values of the Nakagami shaping factor m under power and contention control schemes, 
respectively, in Figs. 8(a) and 8(b). The initial setup in this example is the same as the one 
used for Figs. 7(a) and 7(b), except that the standard variance is a n = 4 dB. When m — 1 
the interfering channel becomes a Rayleigh channel, which is dominated by the log-normal 
shadowing. Whereas, when m = 100 the fluctuations of the channel are reduced significantly 
compared to the Rayleigh fading channel. One fact observed in Fig. 8 is that the interference 
distributions have larger variance and heavier tails when shadow fading is incorporated for both 
control schemes. Interestingly, fading tends to make the interference distribution more heavy- 
tailed than shadowing, i.e., the interference under shadowing has better outage property than that 
under fading. Moreover, the shadow fading has the similar effect for both control schemes. 

VII. Conclusions 

Interference at a primary receiver caused by CR transmitters with power control, contention 
control, and hybrid power/contention control schemes has been characterized. The PDFs of inter- 
ference in the first two cases have been evaluated analytically while, the interference distribution 
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under the hybrid power/contention control has been studied numerically. It has been found that the 
proposed power control and contention control schemes are two effective approaches to alleviate 
interference caused by CR transmitters. The hybrid control scheme causes higher interference to 
a primary receiver, but leads to larger CR coverage as compared to either power or contention 
control schemes. Then, the interference distributions for power and contention control schemes 
have been approximated by log-normal distributions using the cumulant-matching approach 
where the interference PDFs have been obtained with reduced complexity. Furthermore, the 
effect of a hidden primary receiver on the perceived interference has also been investigated for 
the primary receiver. Numerical studies have demonstrated the impact of some CR deployment 
parameters on the resulting aggregate interference under power and contention control schemes. 
It has been shown that increasing the IR radius is an effective way to reduce the interference. 
Moreover, the power control scheme is more sensitive to the IR radius than the contention 
control counterpart. Finally, the impact of shadow fading on the aggregate interference has been 
analyzed as well. 
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Appendix 

A. Derivation of (fl"5T) 

Substituting © and © into CHI), we have 
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</> Y (w)=exp<A7r / f h (h) / / cc (r) 



H 



expiXirj f h (h) J f cc {r) 



R 2 



g(R) 



dr dh 



" (5 (t)) e r P wc dt 



+Att I f h (h) I / cc (r) 



pwc «/ 

rs(«) 



drdh 



iutP m „~h 



dt 



drdh} . (32) 



/f (l - e l ^ p — h ) + iuP max h J (g- 1 (t)) V 
Using ([T]) and (l32l) . the characteristic function (fT5l) is obtained. 

5. Derivation of (TT8T ) 

iC= / / h (/i) / f p { P )^kpd P dh 

pwc / *» \ If /"QO ^ \ 

27rAre~ A7rr dr , (33) 



V -Pmax / fb(h)Vh dh 



H 



2-nr\e 



- A7rr 



where the first equality of (1331 ) holds according to ll2TTl . (TT81) is obatined immediately from (|33l . 
C. Derivation of ([191 ) 

Following similar steps as in [14J, the characteristic function of the aggregate interference can 
be expressed as 



Z— >oo 



(34) 



where 



Q = E U»p B <y)H 
fl 



^iuPg(r)h 



2r 



l 2 -R 2 



A(h) / E 
h Jr. 

f h (h) [ \(1 - q Blh ) + q mh e^(r)h 
H JR 1 



dr dh 



2r 



l 2 -R 2 



dr dh 



1 ~ fc„h + <7mh ^ A(M J e lUJP9{r)h j^Jp 



dr dh. 



(35) 



The integral in the last equality of (1351) can be written as 



lim / f h (h) I e 
1 ^°°Jh Jr 



1 2r 

iu,pg(r)h dr dh = l 



P-R 2 l 2 -R 2 ., u 



f h (h)T(ujph)dh (36) 



where T(uph) is given in flU). Substituting ([35]) and ([36]) into CP}, we obtain (fl9l 



D. Derivation of (|25l) 
From <n~5l) and (124)) . we have 
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K~ I Mh) I /cc(r) 



H JO 



r2 f ir a P max g(R)h \ n nur 



ap u) n rg{R) „ 

max > 1 e^-idt 



pwc 



+ £J B Mh) ! I An 
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The first equality of 071 ) is obtained based on the following fact 



" Qn - 


e aL) = 


" Qn - 


_du n _ 


w=0 


_du n _ 



oj=0 i=Q 



[auj 



In the last equality of (|37l) . the first integral can be expressed as (30 

f h (h)h n dh = e n ^^ 



(37) 



(38) 



(39) 



H 



with ii and a 2 given in © and ©, respectively. Also, the sum of the last two integrals in (|37l ) 
can be simplified as 



f cc {r)r % 



' pwc 



dr + I f cc {r)dr 



na(na — 2) ■ ■ ■ 2 



r pwc « a (27rA) 



Sr( i -' w )-i: 



na((ia — 2) • ■ ■ (no — 2i + 2) 



(27rAr pwc 2 ) 1 



na-2i — A7rr p 
' pwc c 



(40) 



Substituting ([39]) and (@0]) into d37j yields (l25l) . 



£. Derivation of (|28l) 

^y(w) = lim exp {AttDj (j3 (e^-^^) -l)} 

= limexp<A7rD; / f h (h) [ f cc (x) [ — [ exp[iu>p pwc (x)g(r cp (r,6))h] ^—dr d9 dx dh-1 
Jh Jo Jo 211 Jr. 



lim exp { A / f h {h) / f cc (x) 



2r 



exp [«wpp W c(a;)g(r C p(r, 6))h] r — rdr df? dx dh 



(41) 



' H JO Jo JR 

with A = l 2 — R 2 . The first equality in (l4TT) is obtained in the same way as (|34|) and (|35l) . 
Equation (|28|) can be obtained immediately from (14TT) . 
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Fig. 1. System model for CR networks coexisting with a primary network (R = 250 m). 
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Fig. 2. Comparison of interference distributions for power, contention and hybrid power/contention control 
schemes (R =100 m, A =3 user/1 4 m 2 , /3 =4, r pwc = 20 m, a = 4, P max = lW,p=lW, d min = 20 m and 

Hiyb = 30 m). 
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Fig. 3. Log-normal approximation for interference distribution under (a) power control (R =100 m, A =3 
user/10 4 m 2 , (3 =4, r pwc = 20 m, a = 4, P max = 1 W, fj, = and a = 4 dB) or (b) contention control (72 =100 
m, A =3 user/10 4 m 2 , p =4, rf min = 20 m, p = 1 W, fj, = and a = 4 dB). 




Fig. 4. Imperfect knowledge of primary receiver location - the primary receiver is hidden from all CR transmitters 

distributed in the shaded region. 
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Fig. 5. Log-normal approximation for interference distribution with a hidden primary receiver under (a) power 
control {R =200 m, A =3 user/10 4 m 2 , /3 =4, r pwc = 20 m, a = 4, P max = 1 W and r p = Q.5R) or 
(b) contention control (R, =200 m, A =3 user/10 4 m 2 , j3 =4, d min = 20 m, p = 1 W and r p = 0.5R). 
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Fig. 6. Impact of hidden primary receiver on interference distribution for CR networks under hybrid 
power/contention control scheme (R =200 m, A =3 user/10 4 m 2 , f3 =4, a = 4, d ro ; n = 20, p — 1 W, r p = 0.5R 

and Thyb = 30 m). 
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Fig. 7. Impact of various CR deployment parameters on the aggregated interference for CR networks with 
(a) power control (R =100 m, A =3 user/10 4 m 2 , (3 =4, r pwc = 20 m, a = 4 and P max = 1 W) or (b) contention 
control (R =100 m, A =3 user/10 4 m 2 , (3 =4, d mln =20 m, and p = 1 W). 




Fig. 8. Impact of shadow fading on the aggregated interference for CR networks with (a) power control (R =100 
m, A =3 user/10 4 m 2 , (3 —4, r pwc = 20 m, a — 4 and P max = 1 W) or (b) contention control (R =100 m, A =3 

user/10 4 m 2 , (3 =4, d min =20 m and p = 1 W). 



